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Abstract. We investigate the possibility to search for metal-poor high column density (<; 10 23 cm" 2 ) clouds at 
high redshift (z) by using gamma ray burst (GRB) afterglows. Such clouds could be related to primeval galaxies 
which may cause a burst of star formation. We show that a large part of hydrogen is in molecular form in such 
a high column density environment. Therefore, hydrogen molecules (Ha) rather than hydrogen atoms should be 
searched for. Then we show that infrared H2 lines are detectable for metal-poor (^ 0.01 solar metallicity) high 
column density (logA^H [cm -2 ] J; 23.5) systems at high-z without suffering dust extinction. Optical properties of 
dust in infrared could also be constrained by observations of high column density systems. Some possible scenarios 
of producing high column density systems are finally discussed in the context of galaxy evolution. 
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1. Introduction 

The intergalactic and interstellar gas is a main compo- 
nent of the baryonic matter, and occupies most of the vol- 
ume in the Universe. Gas with various column densities 
has been detected as line absorptions of the background 
Quasi-Stellar Objects (QSOs). The strong Lycv absorption 
feature is often used to identify absorbing clouds, which 
are called Lya clouds. In particular, if the Hi column den- 
sity is larger than 2 x 10 20 cm~ 2 , the absorbing system is 
called damped Lya clouds (DLAs) because the continuum 
of background QSOs is completely damped by the Lya ab- 
sorption (e.g., Wolfe et al.^986). 

So far the hydrogen column density has been sampled 
up to log-ZVn ^ 22 (in this paper, we show all the col- 
umn densities in units of cm -2 ). It is difficult to detect 
systems with higher column density, since dust severely 
extinguishes the background objects. Fall et al. (1989) 
show that QSOs behind DLAs tend to be redder than 
those without DLAs. This result strongly suggests that 
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the dust contained in high column density systems extin- 
guishes the background light (but see Murphy & Liske 
2004). If log-ZVu 22, the effects of extinction become 
significant and a large fraction of absorption systems can 
be missed (e.g., Vladilo & Peroux 2005). 

Schaye ( 2001) proposes another possibility of missing- 
high column density systems. He shows that a large frac- 
tion of hydrogen is in molecular form if logiVn <^ 22. Such 
molecular clouds may not be sampled by Lya lines of H I, 
and thus, it is necessary to search for molecular hydrogen 
(Ha). Shibai et al. 1|20011 hereafter S01) investigate the 
observability of H2 in absorption for high column density 
systems. They focus on near-infrared (NIR) vib-rotational 
lines, because dust extinction is less severe in the NIR 
than in the optical and ultraviolet (UV). Considering typ- 
ical QSOs as background sources, they show that the 
high column density systems with an H2 column density 
of log-ZVn 2 <; 23 could be observed without being extin- 
guished by dust if the dust-to-gas ratio is less than 1 per- 
cent of the Galactic value. 

However, QSOs are usually sampled in the optical 
wavelength, where dust extinction is severe. Ellison et 
al. (2001) investigate a sample of radio-selected QSOs 
and conclude that dust-induced bias in optical samples 
is small. Therefore, the significance of extinction bias is 
still debated. In this paper, we consider the possibility 
to use gamma-ray burst (GRB) afterglows as background 
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sources. Since GRB samples are defined in wavelength not 
affected by dust extinction, we expect less extinction bias 
in GRB samples than in QSO samples. Indeed some DLAs 
are found to lie in the sightline of GRBs (e.g., Vreeswijk et 
al. I2004j) and some GRBs are surrounded by high column 
density gas (e.g., Piro et al. I2005[) . Therefore, it may be 
possible to use NIR continua of GRB afterglows to detect 
NIR H 2 lines. 

Recently, Inoue et al. (2006) have investigated the pos- 
sibility of detecting the absorption lines originating from 
protostellar cores by using background GRB afterglows. 
They find that radio observations of some metal absorp- 
tion lines of GRB afterglows are useful to detect high- 
z protostellar gas. However, it is still interesting to di- 
rectly detect H2 because H2 is considered to play a central 
role in star formation in metal-poor environments (e.g., 
Mizusawa et al. I2004f) . Therefore, in spite of its difficulty, 
the possibility to detect IR H 2 lines is worth investigating. 

Since H2 vib-rotational lines are generally weak, large 
infrared (IR) telescopes are necessary. The strongest vib- 
rotational H2 lines are normally those around the wave- 
length of A ~ 2 ^m. Those lines shift to the mid-lR (MIR) 
regime if we observe high-z clouds. In MIR, it is essen- 
tial to avoid the emission of atmosphere and telescope. 
Thus, a cooled space telescope is ideal. There are indeed 
suitable future missions planned currently such as the 
Space Infrared Telescope for Cosmology and Astrophysics 
(SPICA) 1 . Herschel 2 is also a planned space mission in 
IR (but not cooled), and the same kind of targets can be 
accessible if the background sources are luminous enough. 
In NIR (A ^ 5 /im), the James Webb Space Telescope 
(JWST) 3 is more sensitive than the SPICA, and is ap- 
propriate for H2 absorption at z 1.5. 

First, in Sect.|2]we explain the method to calculate H2 
absorption line luminosity. Our results are shown in Sect. 
13 and some relevant physical and observational issues are 
discussed in Sect.0] Our conclusions are given in Sect.[SJ 

2. H2 absorption lines 

2.1. Absorption line flux 

We follow SOI for the formulation concerning IR H 2 lines. 
We assume a uniform, cool gas cloud satisfying k-QT cyi <C 
hv, where fee is the Boltzmann constant, T cx is the exci- 
tation temperature, h is the Planck constant, and v is the 
frequency of an H2 line. 

The optical depth of the line absorption, rii nc , is esti- 
mated as 

riine ^(^W^_L, (I) 
8tt \g e J Av 

where subscripts u and t indicate the upper and lower lev- 
els of a transition, g u and gi are the degeneracy of each 
state, respectively, A u i is Einstein's A coefficient, Ni is 

1 http:/ /www. ir.isas.ac.jp/SPICA/index. html 

http:/ /astro. estec.esa.nl/SA-general/Projects/Herschel/ 
3 http://ngst.gsfc.nasa.gov/ 



the column density of the molecules in the lower state, 
and Av is the line width in units of velocity. If the instru- 
mental wavelength resolution is larger than the physical 
line width, we should take Av = c/R, where R is the 
resolution defined by R = A/AiA (AiA is the wavelength 
resolved by the instrument). 

Next, we consider extinction by dust grains, whose op- 
tical depth at the wavelength A is denoted as 7d U st(A). 
This is estimated as (Hirashita et al. 2 0050 

Tdust(A) = ^TO H A r H^ , (o'd(A)/m d ) , (2) 

where pi is the mean atomic mass per H atom (assumed 
to be 1.4), rein is the mass of a hydrogen atom, ATr - is the 
column density of hydrogen nuclei, T> is the dust-to-gas 
mass ratio, and (<7d(A)/md) is the cross section of dust 
per unit dust mass as a function of wavelength. 

Finally, we obtain the absorption line flux with extinc- 
tion, igj£, as 

= S v Av[l - cxp(-T linc )] cxp(-T dust ) , (3) 

where Av (= v Av/c) is the line width in units of fre- 
quency and S v is the continuum flux (per frequency) of 
the IR source behind the cloud. 

We assume that all the molecular hydrogens are in the 
levels J = 0, I of the vibrational ground state (v = 0), 
since we are interested in the clouds whose excitation 
temperature is much lower than the excitation energy. 
The lines are summarized in SOI. The ortho vs. para 
ratio is assumed to be 3 : I; i.e., N\ = (3/4)7Vh 2 and 
Ao = (1 / A) Nb_ 2 , where Nj is the column density of H2 in 
the level J, and Nn 2 is the total hydrogen column density. 

Observationally it could be convenient to present the 
equivalent width, since it is not sensitive to the line width 
(or the spectral resolution). The equivalent width W\ 
given in wavelength units is expressed as 

A 2 

W x = — [l-exp(--n inc )]Ai/. (4) 
c 

We calculate W\ as a function of ATh 2 as shown in Fig. 
^ We should note that not only the equivalent width but 
also the suppression of continuum level by dust extinction 
is important for the absorption flux (Eq. |3j) . 

2.2. Extinction models 

The optical depth of dust, Td us t, is proportional to the 
dust-to-gas ratio, and the hydrogen column density, if we 
fix (<7d(A)/rad) as a function of A (Eq. |2J), which is taken 
from Hirashita et al. l|2005j) . We adopt the Galactic extinc- 
tion curve, which is consistent with Mathis (jTSHQJ, unless 
otherwise stated. A different extinction curve is examined 
in Sect. 14.21 We assume that the Galactic dust-to-gas ra- 
tio, 2?q, is 6 x I0~ 3 , since we should adopt the values 
in Hirashita et al. I|2005[l for consistency. We define the 
normalized dust-to-gas ratio, 

K = V/V Q . (5) 
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Fig. 1. Equivalent width (EW) for various H2 lines as a 
function of H2 column density. The solid, dotted, dashed, 
and dot-dashed lines represent the line wavelengths of A = 
2.12, 2.22, 17, and 28 fim, respectively. 



2.3. Background sources 

We focus on GRB afterglows as one of the most lumi- 
nous sources in the Universe. If some dense clouds lie in 
the line of sight of GRBs, we can have a chance to de- 
tect these foreground clouds in absorption. For the flux of 
GRB afterglows, we adopt a commonly used model with 
synchrotron radiation from a relativistic shock (Sari et al. 
1998] Ciardi & Loeb l2l)l)Tj|) . Here we apply a set of values 
adopted by Inoue et al. I|2004(l : a magnetic energy frac- 
tion of 6b = 0.1, an electron energy fraction of e c = 0.2, 
a spherical shock energy of E — 10 52 ergs, an ambient 
gas number density of n = 10 cm , and a power-law in- 
dex of the electron energy distribution of p = 2.5. For the 
timescale, we put t = 1/24 day (1 hour), and the obser- 
vational detection limits are also calculated with 1-hour 
integration. 

We do not include the redshift dependence of GRB 
luminosity. GRBs may be brighter at high z than the local 
Universe (e.g., Lloyd-Ronning et al. 120021 Yonetoku et al. 
2004). In this case, detection of high-z absorption lines 
becomes easier. 



3. Results 

We consider a velocity width of Av = 100 km s _1 , moti- 
vated by the future cooled space telescopes such as SPICA 
with R ~ 3000. Accordingly, the frequency resolution be- 
comes = v/R. Then, S u Ai/ is estimated as 

s ^ = 1 .„ x !„-..(_!_) (_A_)-'(_^)- 1 

WnT 2 . (6) 



We assume that the all the hydrogen atoms are in the 
molecular form (i.e., = 2Ah 2 ). This assumption is 
discussed in Sect. 14.11 

We note that S01 scale the extinction by Z (metallicity 
normalized to the solar value). We have confirmed that if 
we adopt k — Z , the difference in Td us t less than a factor 
of 2. As shown by S01, the 2.12 /mi absorption line is the 
strongest, and the dust obscuration does not affect the 
line detection if k < 0.01. Therefore, we concentrate on 
the 2.12 fim line with k = 0.01 in this section. 

As stated by S01, it is quite difficult to detect absorp- 
tion lines where optical depth is less than 0.01, since an 
extremely high signal-to-noise ratio would be required to 
detect an absorption with rii nc < 0.01. The H2 column 
density which satisfies rii ne = 0.01 is listed in S01 for vari- 
ous lines. For the 2.12 /zm line, Nb_ 2 should be larger than 
4 x 10 23 cm -2 to satisfy rii nc > 0.01. Therefore, if we use 
the NIR H2 lines, we should target the absorption sys- 
tems with log7VH 2 ^ 23.5. However, we should note that 
this condition is dependent on the line width. The sound 
speed of H2-rich cold clouds is generally much less than 
100 km s . Therefore, if we had a facility sensitive enough 
with a resolution large enough, we could observe an ab- 
sorption system whose H2 column density is smaller than 
1Q 23.5 cm -2 

In Fig. [3 we show the 2.12 /zm absorption fluxes as a 
function of -/Vh 2 for various background fluxes. The typi- 
cal detection limit expected for future IR cooled telescopes 
such as SPICA is also shown (5 x 10~ 21 W m~ 2 for 5 a 
detection with an integration of ~ 1 hour; Ueno et al. 
2000] S01) as a representative observability of future fa- 
cilities. With future space missions, high column density 
systems with log 7Vh 2 ~ 23.5-25 can be investigated if the 
background source luminosity is <; 10 mJy. 

In order to show the flux level of GRB afterglows, 
we present Fig. |3] The solid and dotted lines show the 
flux (S„) at the restframe wavelengths of A = 2.1 /im 
and 17 /zm (i.e., the observational wavelengths of 2.1(1 + 
^grb) M m an d 17(1 + zgrb) A™), respectively, where zgrb 
is the redshift of the GRB afterglow. We observe that the 
afterglows are more luminous at 2.1(1 + zgrb) A™ than 
at 17 (1 + zqrb) M m - In reality, we should observe at 
2.1(1 + z abs ) A*m or 17(1 + z abs ) /zm, where z abs is the 
redshift of a target absorption line system, but if z a b s has 
the same order of magnitude as 2grb , the absorption flux 
does not change very much (the dependence on other pa- 
rameters assumed in Section 12.31 is larger) . Considering 
that the 2.1 /zm absorption line is stronger than the other 
lines, observations at 2.1(1 + z a bs) A 4111 is favorable to trace 
high column density (logiVH 2 ~ 23.5-25) systems (see S01 
for the other lines). 

The dependence on the dust-to-gas ratio is also shown 
for the background flux of 10 mJy (Fig.0J. The solid, dot- 
ted, dashed, and dot-dashed lines represent the absorption 
line fluxes with k = 0.0001, 0.001, 0.01, and 1, respectively. 
We see that if k > 0.1, we can detect the NIR H2 line if 
logA^H 2 23. 
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Fig. 2. 2.12 /xm line absorption flux as a function of the 
column density of hydrogen molecules. Various continuum 
fluxes of the background sources are examined. The solid, 
dotted, dashed, and dot-dashed curves correspond to the 
continuum fluxes of 30, 10, 3, and 1 mJy, respectively. 
The horizontal thin solid line shows the expected detection 
limit of future cooled space telescopes such as SPICA (5 
a for an integration of ~ 1 hour). 
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Fig. 3. Observed flux of a gamma-ray burst afterglow as a 
function of redshift. The solid and dotted lines correspond 
to the restframe wavelengths of 2.12 and 17 /zm, i.e., ob- 
servational wavelengths of 2.12(1 + z) and 17(1 + z) /xm, 
respectively. 

4. Discussions and observational strategies 

4.1. High column density systems in cosmological 
scenarios 

4.1.1. Protogalactic disks 

Collapsed gas in a dark matter halo may form a disk 
because of the initial angular momentum. If we observe 



Fig. 4. Same as Fig. |2 but for various dust-to-gas ra- 
tio (k is the dust-to-gas ratio normalized to the Galactic 
value). The continuum fluxes of the background sources is 
assumed to be 10 mJy. The solid, dotted, dashed, and dot- 
dashed lines correspond to k = 0.01, 0.1, 1, and 0.001, re- 
spectively. The lower dotted line represents the result with 
the H05b extinction curve, while the other lines adopt the 
Galactic extinction curve. 

a GRB in a nearly edge-on disk, we could identify the 
disk as a high column density system. If a galactic disk 
has a typical circular velocity of 100 km s _1 , the typi- 
cal mass (Mdisk) and radius (rdisk) of the disk should be 
Mdisk ~2x 10 9 Mq and rdisk ~ 6 kpc if it forms at z ~ 3 
(Hirashita ct al. 2003 ) . With these mass and size, the typ- 
ical hydrogen column density along the galactic disk is 
~ ?Wdisk- Using the relation Mdisk ~ 7rr disk- ff / xm H^H, 
where H is the disk thickness, we obtain 

1 ir>23 ( M disk \ ( 7"disk \ 1 

" Hrdisk ~ 1Xl ° Ux10»mJ U]^J 

x (rof^) lcm "- < 7 > 

Therefore, a typical galactic disk could be observed as a 
high column density system if a GRB lies within or be- 
hind the disk and if a significant fraction of the gas is in 
molecular form. 

As stated in Sect. it is required that the dust- 
to-gas ratio is Si 10~ 2 times the Galactic value (i.e., 
T> 6 x 10~ 5 ). With the typical disk mass mentioned 
above (2 x 10 9 M Q ), the corresponding dust mass is 
Miust ^ 10 5 Mq. We estimate the timescale of the dust 
enrichment up to this level. Galactic disks observed as 
DLAs can be considered as objects with a star formation 
rate of - 0.1 - 1 M yr" 1 (Hirashita & Ferrara 12711751 
hereafter HF05). Using the Salpeter initial mass function, 
we obtain the supernova rate of ~ 10~ 3 — 10~ 2 yr -1 . 
Assuming that around 0.5 M Q of dust is formed in a SN 
(Todini & Ferrara EUJUl Nozawa et al. I20«3|) . we obtain a 
typical dust formation rate of ~ 5 x 10~ 4 — 5 x 10~ 3 M Q 
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yr -1 . Thus, in order to accumulate the dust mass esti- 
mated above, a timescale of 2 x 10 7 — 2 x 10 8 yr is re- 
quired, which is much larger than the lifetime of the GRB 
progenitors (Heger et al. G2003 ) . 

It is observationally known that some DLAs have dust- 
to-gas ratio of k ~ (Pettini et al. 119971 Ledoux et 
al.|2D03. HF05). This observational evidence strongly sug- 
gests that galactic disks with n ~ 10~ 2 exist at high red- 
shift. 

4.1.2. Supergiant molecular complexes 

In some starburst galaxies such as the Antennae (Wilson 
et al. I2TM)> and Arp 229 (Casoli et al. 1999), super- 
giant molecular complexes are observed. In those galax- 
ies, molecular gas of ~ 10 9 Mq is concentrated in a 1- 
kpc 2 area. This indicates an H2 column density as large as 
~ 10 23 cm -2 (and a hydrogen number density of nn ~ 30 
cm -3 ). Arp 220 seems to have a comparable amount of gas 
in more concentrated region (Sakamoto et al UMty . Thus, 
the H2 column density in a starbursting region can exceed 
several times 10 23 cm~ 2 . This column density meets the 
requirement to detect H2 in absorption. 

HF05 have recently developed a method to calculate 
a probable range of /h 2 under a given range of physical 
quantities concerning ISM. Their assumption is that /h 2 is 
determined by the equilibrium between formation on dust 
grains and destruction by dissociating photons. The prob- 
ability density function of log/jj 2 , p(log/n 2 ) is calculated 
by using the method described in Appendix B of HF05. As 
a first guess, we start from the ranges of quantities derived 
for H 2 -detected DLAs by HF05: 1.5 < logn [cm" 3 ] < 2.5, 
1.5 < logT [K] < 3, and 0.5 < logx < 1-5, where n is the 
number density of gas particles, T is the gas temperature, 
and x is t ne UV interstellar radiation field intensity nor- 
malized to the Galactic value taken from Habing (1968). 
We adopt k = 0.01 for the dust-to-gas ratio (Sect.[3J). In 
Fig. we show p(log/n 2 ) for various Ah- From Fig. CD 
we observe that the gas becomes almost fully molecular 
for log Ah 23 even if the dust-to-gas ratio is only 1 
percent of the Galactic value (k ~ 0.01). This is consis- 
tent with the conclusion by Schaye (200T)l that the large 
part of hydrogen is in molecular form at log Ah <; 22. 
We stress that at k ~ 0.01 H2 formation is accelerated 
and possibly leads to a burst of star formation (Hirashita 
& Ferrara 2002). Therefore, our proposed observation has 
an importance that we trace the reservoir of gas clouds 
which might cause a starburst. Indeed Iliev et al. (2006) 
investigate the fate of those clouds proposed by HF05 and 
conclude that such clouds can finally collapse even in the 
presence of UV radiation. 

4.2. Extinction curves 

Some observations show that the extinction curves of high- 
z galaxies are different from local ones (Maiolino et al. 
2004). Therefore, it is interesting to examine the extinc- 
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Fig. 5. Probability distribution function of the molecular 
fraction (log/n 2 ) for various hydrogen column densities. 
The solid, dotted, dashed, dot-dashed, and dot-dot-dot- 
dashed lines correspond to log Ah (cm -2 ) = 21, 22, 22.5, 
23, and 23.5, respectively. The dust-to-gas ratio normal- 
ized to the Galactic value is assumed to be k = 0.01. 

tion curve derived for high-z galaxies. We adopt the the- 
oretical curve of Hirashita et al. I|2005)) . We adopt Model 
b (based on the dust production in an unmixed SNe II; 
Nozawa et al. l2TT0~gj) of Hirashita et al. I)2005|) . which nearly 
reproduce the observational high-z extinction curve of 
Maiolino et al. 1|2004[1 • This model is labeled as H05b (see 
Takeuchi et al. 120051 for the NIR part of the extinction 
curve). 

In Fig. 0] we show the absorption flux of the 2.12 /an 
line for the H05b extinction curve by the lower (thin- 
ner) dotted line, where the dust-to-gas ratio is assumed 
to be k — 0.01 (0.01 times the Galactic value). The upper 
(thicker) dotted line corresponds to the Galactic extinc- 
tion curve with the same dust-to-gas ratio. The difference 
in the extinction curve does not affect the absorption line 
flux for log Ah 2 ^ 23, since the optical depth of dust is 
small, but the difference in the extinction becomes sig- 
nificant for logAn 2 ^ 23.5. In H05b the system becomes 
more opaque for logVn 2 ^ 24 than in the Galactic ex- 
tinction case. This sensitive dependence on the extinction 
could also be used to constrain the IR extinction curves 
at high z (see also Vergani et al. 2004 ). 

4.3. Number Estimate 

There is no observational constraint on the number of high 
column density systems with log A^h <^ 22. One simple way 
to estimate the number is to extrapolate the distribution 
of Ah sampled up to log Ah ~ 22, although Vladilo & 
Peroux (2005| argue that the distribution should not be 
extrapolated in the range of column densities not sampled 
by the observation. After correcting the distribution for 
dust extinction, Vladilo & Peroux l|2005|l indicate that 
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log/(7v" H ~ 10 22 cm- 2 ) 23.5, where f{N H )dN H dX 

represents the number of absorbers along a random line 
of sight with hydrogen column densities between Nn and 
TVn+diVH and redshift paths between X and X+dX (e.g., 
Peroux et al. l20TI5|) . By assuming /(iV H ) oc N H P with = 

1.5 (Vladilo & Peroux l2nn5| l. we obtain log/(JV H ) 25.8 

at log7V H ~ 23.5. At z = 2, X ~ 4.4 and if we take the 
interval of Ns as ~ 10 23 cm -2 , we obtain the number of 
absorbers in a random line of sight as ~ 10 -2 8 . Therefore, 
around 600 lines of sight are necessary to detect a high 
column density object. If a large part of those sources are 
favorable for H2 formation as suggested in Sect. 14.11 we 
can detect H2 lines every ~ 600 GRBs. 

For example, the sample of Schmidt (2001) consists of 
1391 GRBs in a period of 5.9 yr. Their analysis indicates 
that about 500 GRBs are around z ~ 2. Therefore, it is 
possible that about 100 GRBs are identified around z ~ 2 
per year. Therefore, we could find 0.2 GRBs associated 
with high column density systems per year. If we consider 
an improvement of the sensitivity in the future, we expect 
a larger number. 

5. Conclusions 

In this paper, we have investigated the possibility to de- 
tect IR H2 absorption lines from dust-poor high column 
density (logiVn ^ 23.5) clouds at high-z. H2 is a unique 
tracer of high column density clouds, because hydrogen 
is considered to be in molecular form even at extremely 
low dust-to-gas ratio (n ~ 0.01, where K is the dust-to-gas 
ratio normalized to the Galactic value) (Sect. 14. If) . We 
have shown that if K < 0.1, we can observe H2 absorption 
lines of gas whose column density is logA^ ii 23.5. In 
particular, if k ^ 0.01, clouds with log7\?H 2 ~ 23.5-25 are 
accessible (Sect.|3 Fig.0J. Since near future missions such 
as SPICA will access GRB afterglows up to z ~ 2 in the 
IR regime, the search for metal-poor (or dust-poor) galax- 
ies at z ^ 2 is an important scientific target of near future 
IR telescopes fSect. l4~3l . We should note that some bright 
GRBs at high-2 can be utilized to trace higher- z primeval 
galaxies. 
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